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Comparing the results of exact quantum calculations and those obtained from the EBK-like 
quantization scheme of Silvestrov et al [Phys. Rev. Lett. 90, 116801 (2003)] we show that the 
spectrum of Andreev billiards of mixed phase space can basically be decomposed into a regular and 
an irregular part, similarly to normal billiards. We provide the first numerical confirmation of the 
validity of this quantization scheme for individual eigenstates and discuss its accuracy. 
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For quantum systems, whose classical analogues have 
mixed dynamics, the separation of the spectra into reg- 
ular and irregular parts goes back to PercivaA Here 
the term 'regular' refers to such energy levels, which 
correspond to quantised invariant tori, while 'irregu- 
lar' refers to those associated with the chaotic part of 
the phase space. Semiclassical methods, such as EBK 
quantization 2 , have given a deep understanding of the 
influence of the classical dynamics on the quantum spec- 
trum (see eg2A£ and references therein). 

Recently a new type of quantum dot system, compris- 
ing normal-superconductor (N-S) interface has attracted 
considerable attention. These systems are commonly 
called " Andreev-billiards" (ABs)££ The name derives 
from the specific scattering process taking place at the N- 
S interface, namely the Andreev-reflection 8 , whereby an 
impinging electron-like quasiparticle with energy e (mea- 
sured from the Fermi energy E-p) is coherently scattered 
into a hole (and vice versa) if e is smaller than the su- 
perconducting gap A (for details see eg 7 ). Unlike spec- 
ular reflection, this scattering process is accompanied by 
the (approximate) reversal of all velocity components, 
thereby giving rise to a peculiar classical dynamics. In 
a recent paper— we have presented evidence that in gen- 
eral the phase space of these systems is mixed. Thus the 
question naturally arises, of whether one can perform a 
similar separation of the spectrum into regular and irreg- 
ular parts as in normal systems. While the influence of 
different regions of the mixed phase space of an isolated 
normal dot on the density of states of the corresponding 
ABs has been adressed beforehand^ to our knowledge no 
study has tried to answer the above question regarding 
the individual energy levels and taking into account the 
peculiar dynamics of the whole N-S system. Our work is 
aimed to be the first step towards the answer by studying 
a simple yet nontrivial 2D example. 

In our study we shall use semiclassical and quan- 
tum mechanical tools to identify the regular eigenstates. 
Quantum mechanically, these systems can be described 
by the Bogoliubov-de Gennes equations^. From semi- 
classical point of view, this implies that the EBK quan- 
tization has to be generalized to the case of spinor wave 
function. Such a multicomponent semiclassical theory 



for the N-S systems has been developed in Ref. 12 and it 
was shown that the EBK quantization of normal systems 
can also be generalized to integrable N-S systems. In- 
deed, in Ref™ a good agreement between semiclassical 
and quantum calculations has been found for an inte- 
grable AB. Another important idea regarding the quan- 
tisation of N-S systems was introduced in RefsJ^A^ for 
cases where for e > the classical dynamics is strictly 
speaking non-integrable, but there is an adiabatic invari- 
ant. At zero magnetic field the invariant is the time T(e) 
between subsequent Andreev reflections. According to 
this approach, for the purpose of semiclassical quantiza- 
tion, one can consider the curves of constant T at e = 0. 

In our work we benefit from both of the studiesi 2 ^. As 
it will be shown below, the dynamics in certain regions of 
the phase space is quasi-integrable. By quasi-integrable 
we mean regions which contain mostly tori, on which the 
dynamics is similar to that of integrable normal systems. 
Therefore we expect that for these islands of regular mo- 
tion the results of Ref™ apply. However, the analyti- 
cal calculations can be performed only in the adiabatic 
approximation^, whose accuracy will also be discussed. 

The model system we used in our calculations is the 
Sinai-Andreev (SA) billiard2ii&. It consists of a Sinai- 
billiard-shaped normal dot and an attached (infinite) su- 
perconducting lead as shown in Fig. [He). Classically, in 
the phase space of this model we have found a large sta- 
bility island (see Fig.UJa)). Its existence can be expected, 
since it is centered on the shortest unstable periodic or- 
bit of the isolated normal dot (which corresponds to the 
motion along the bottom wall) and the presence of the 
superconductor, i.e., Andreev-reflection can stabilise this 
orbit. To prove this, we consider the dynamics on the 
Poincare section (PS) which we define in the following 
way: we record the position y and the tangential veloc- 
ity component v y in units of v e = VfV^ + £ when the 
Andreev reflection results in a departing electron quasi- 
particle, ie whenever a hole impinges on the N-S inter- 
face (here e = e/Ep, vp is the Fermi velocity and y is 
measured from the lower edge of the interface). Suppose 
now that an electron departs from the interface and after 
returning to it becomes Andreev-reflected. Using sim- 
ple geometrical considerations and taking into account 
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the Andreev reflection lawi£ii£ one finds that for y <C 1, 
v y <C v e the linearised equations of motion for the phase 
space coordinates y, v y = v y /v e of the quasiparticle is 
given by the stability matrix 

M (e)-( 1 - 2d/R 2d (l-d/R)\ a) 
M(6) -{ 1 2/R 1 (2d/R-l) ) (L) 

where 7 = \/l + ej \J\ — e and the geometrical param- 
eters <i, R are defined in Fig. HJc). The motion of the 
emerging hole can then be described by M(—e). How- 
ever, by construction the Poincare map consists only of 
the starting coordinates of electron trajectories, there- 
fore the motion in the (y,v y ) plane is given by M e h = 
M(—e) M(e). The periodic orbit is stable, if the trace of 
the matrix M e h is less than 2, which gives the condition 
< d < R/(j — 1). Since e <C Ep in our calculations, 
1/(7— 1) ~ 1/e ^> 1 and the periodic orbit is stable for 
wide range of values of the parameters d and R. 

Part of the PS around the stable periodic orbit is shown 
in Fig.^a). As expected, we find numerically that there 
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FIG. 1: (colour online) Enlargement of a part of the PS 
around v y = (a) (for the geometrical and other parameters 
see—). Each dot represents a starting point of an electron 
trajectory. The red curves show the ellipse-approximation 
given by Eq. Q, for different T s. The blue curve denotes an 
another, thin regular island. An alternative view of PS show- 
ing bands of approximately constant T e h denoted by different 
colours (b). Only bands with T e h < 6a /vf are shown. One 
of the tori shown in (a) is projected onto the real space (c). 
The electron trajectories are blue, the hole ones are red. 

are truly invariant curves around the stable periodic orbit 
forming an island of regular motion. This regular island, 
which can be clearly seen in the middle of Fig.^Ja) is sur- 
rounded by intermittent- like regions. Beyond there are 
chaotic seas intervowen with other intermittent regions 
and there is at least one another smaller regular island. 
As pointed out in our previous paper—, this structure is 
due to the interplay of non-exact retracing for e > and 
the presence of the so called "critical points" which sep- 
arate the normal and the superconducting segments of 
the billiard boundary. 

Because of its importance in the semiclassical quanti- 
zation, we now focus on the energy dependence of the 



time T e h(y,v y ,e) of the electron- hole orbit between two 
section with the PS in the regular island. Following the 
electron-hole orbits which trace out the invariant curves 
Ci nv in Fig. nja), we have found numerically that (a) 
the variation of T e h(y, v y ,s) corresponding to subsequent 
points on the invariant curves is only 0(e 3 ), ie T e h is an 
adiabatic invariant (b) T e h = 2T(y,v y ,e = 0) + 0{e 2 \ 
where we denoted by T the time between subsequent 
Andreev-reflections for e = 0. An important consequence 
of the above observations is the following: if we denote 
by Ct the T(y,v y ,e = 0) = To curves on the PS, then 
considering any point (y,v y ) on Ct , for finite e the invari- 
ant curve Ci nv which contains this particular point will 
always be in 0(e 2 ) vicinity of Ct - This will be important 
later on when we discuss the semiclassical quantization of 
this regular island: it can be showni2ii4, that one of the 
action integrals to be calculated equals the area enclosed 
by Ci nv , which can be then approximated by the area en- 
closed by Ct , since the difference is of higher order than 
e <C 1. A similar conclusion has been drawn in Ref™. 

An alternative view of the PS can be obtained by de- 
noting with different colours those regions, which cor- 
respond to approximately constant T e h. These regions, 
which usually appear as narrow bands in Fig. W[b) are 
often not easily recognizable on the PS since they do not 
always show ordered pattern such as those indicating the 
presence of intermittent-like motion around the quasi- 
integrable island in Fig. [Ha). (The large dark region in 
the middle of Fig. Gib) corresponds to the regular island 
and the surrounding intermittent region in Fig. IHa).) 
As we will briefly discuss later, the importance of these 
bands lies in the fact that some of them can be associated 
to certain eigenstates. 

Before turning to the adiabatic quantization, we note 
that for a different W, and larger value of e, in the regu- 
lar island we have also found structures resembling very 
much the secondary island chains known from normal 
KAM systems. It would be an interesting future project 
to investigate the properties of the classical dynamics in 
more detail to explore the exact nature of the similari- 
ties with KAM systems. However, for the parameters^ 
used in the present study, these "secondary islands", if 
they exist, are much smaller than Planck's constant h 
and thus do not play any role in our further discussion. 

We now proceed with the adiabatic quantisation^ of 
the regular island shown in Fig^a). The semiclassical 
energies can be obtained by quantising the action inte- 
grals Ii = l/2it § c . pdq, i = 1,2, calculated along the 
Ci irreducible closed contours on the adiabatic tori. The 
curve Ct can be chosen as the integration contour C\ 
and therefore I\ equals the enclosed area on the PS. One 
can show that for y, v y <C 1 the curve Ct is semi-ellipse 
given by the following equation: 

(L-2d)R 2 = Fy 2 + Gv y y + Hv 2 y (2) 

where F = (2d + R), G = 4d(d + R), H = FG/4 and L = 
v-pTo. Comparison of the above curve for different Xbs 
and the numerically-calculated PS of the corresponding 
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tori for e = 0.0105 is shown in Fig.QJa). One can observe 
that the agreement is very good for the inner tori, while 
for those closer to the border of the regular island there 
is a small but noticeable deviation since the assumption 
^{iyCl does not hold. Quantising the area enclosed by 
Ct one finds the condition 



1 (L - 2d)R 
8 y/d(d + R) 



Pf = h [ m 



(3) 



where m = 0, 1, 2, . . . and the Maslov index is 3/4 since 
one caustic and one hard wall is encountered along the 
integration contour. 

The second action integral I2 in the adiabatic approxi- 
mation is to be calculated along a self-retracing electron- 
hole orbit. Since the particles move ballistically inside 
the normal dot, the quantisation condition reads: 

(p e — Ph)L = 27rh[n + — arccos(s/A)]. (4) 
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Here p e (e) = £>f\A + £ and ph(s) = Pe(— s) are the mag- 
nitudes of the electron's and hole's momentum respec- 
tively, n = 0, 1, 2, . . ., and with the arccos(£/A) term we 
take into account the phase shift due to the Andreev- 
reflection. Combining Eqs. (JjJJ and (J3J) one arrives at an 
implicit equation for the eigenvalues e nm which can be 
solved numerically for each value of n and m. 

Before comparing the semiclassical and exact quantum 
calculations, we briefly discuss the accuracy of the adi- 
abatic approximation. There are two sources of error: 
the first is in the derivation of Eqs. © and (|2J) where 
we considered the invariant surfaces in the phase space 
for e = 0. However, due to the non-exact retracing, the 
dynamics is different for finite e. Since in the quantum 
system the excitation energies e nm are always finite, one 
should calculate the action integral I\ using the invariant 
curves Ci nv , whereas in case of I2 one should take into ac- 
count that for i > the hole does not retrace exactly the 
path of the electron, and one should choose the integra- 
tion contour accordingly. It turns out however, that for 
this regular island the results for I\ and I2 do not change 
in first order of e even if we take into account the effects 
of finite e on classical dynamics. Namely, as we have 
already pointed out, the difference between the areas en- 
closed by a Ci nv and a corresponding Ct is only of 0{e 2 ) 
and therefore I\ is accurate to first order in e. Consider- 
ing I 2 , we checked the accuracy of Eq. (@J) by taking as an 
integration contour C2 a non-retracing electron-hole orbit 
and then we closed the contour on the PS. After lengthy 
calculations we have found that the first correction to the 
result given by Eq. Q is of the order of e 3 . 

The second source of error is that in order to calculate 
1\ we approximated the curves Ct by the semi-ellipses 
given by Eq. This certainly introduces inaccuracy in 
case of adiabatic tori close to the border of the regular is- 
land, for which the conditions |/,^ <C 1 do not rigorously 
hold. Therefore we checked our analytical results in the 
following way: first we numerically determined those Ct 



curves on PS for which the enclosed area is Ti(m + 3/4), 
m = 0, 1, .... Then reading off the Tq values correspond- 
ing to these curves and using Eq. (@J) we re-calculated 
the semiclasssical energies. The results are summarized 
in Fig.|2ta),(b) and will be discussed below. 

We now compare the results of the quantum mechani- 
cal and semiclassical calculations and show that the semi- 
classical energies agree remarkably well with the exact 
quantum ones. The quantum treatment of the system 
is based on the Bogoliubov-de Gennes equations^ and 
is briefly described inS. It is assumed that the super- 
conducting pair potential A is constant inside the lead 
and zero in the N region! and we work in the regime 
Sn <C Et « A« E'pi where S n is the mean level spac- 
ing of the isolated normal dot and Et is the Thouless 
energ3*i£. The area of the regular island on the PS is ~ Ah 
and therefore we expected 8 regular eigenstates in the 
spectrum, corresponding to quantum numbers n = 0,1, 
m = 0,1,2,3 (For the parameters^ used in our calcu- 
lation Eqs. @, (0 do not have real solution for larger 
n, while the values of m are limited by the size of the 
stability island). 
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FIG. 2: (colour online) Quantum eigenenergies (green 
squares), the results of the adiabatic quantization using 
Eq. © (red stars) and by taking into account the exact shape 
of Ct (blue triangles) for n = (a) and n = 1 (b). The arrows 
show which quantum and semiclassical energies correspond to 
each other. The bars show the magnitude of the mean level 
spacing £ns = Sn/2. The square modulus |^(r)| 2 of the elec- 
tron components of the wave functions corresponding to cer- 
tain eigenenergies as indicated by the letters c, d in Fig. EI a) 
are seen in (c),(d). The Kr)| 2 for the state at e/A = 0.3202 
corresponding to a band of approximately constant T e h start- 
ing at y = 0,v y = ±0.4 in Fig. ^b) is shown in (e). 

Comparison of the semiclassical predictions for n = 0, 
m = 1 ... 4 and the exact quantum mechanical eigenval- 
ues lying roughly in the same energy range can be seen 
in Fig. Efa) and (b). This shows that there are twice as 
many quantum eigenenergies as semiclassical ones and at 
first sight it is not obvious which of the eigenstates should 
be considered as regular ones. To identify the eigenstates 
which correspond to quantized tori we now examine the 
eigenfunctions. Similarly to normal billiards^ we have 
found that the wave function of certain eigenstates show 



strong localization (both the electron and the hole com- 
ponents) onto classical objects (see alsc&i^). This can 
be either a torus (compare Fig. 0(c) and Fig.^c)), or a 
bunch of trajectories (see eg Fig. 0(e)) corresponding to 
a band of approximately constant T e h in Fig.^b). 

Computing eg the eigenfunctions belonging to the 
eigenvalues shown in Fig. 0^a) reveals that only four of 
them are localized (as in Fig. 0^c)), three others display 
an apparently random interference pattern and cover the 
normal dot in roughly uniform way (as in Fig. 0d)), 
while one of them is of intermediate nature. (The pic- 
ture is very similar in case of Fig. 0b)). As an example, 
we show the electron component of those two eigenstates 
energies of which are very close to the n = 0, m = semi- 
classical one (see Figs. 0c), (d)). One can clearly see that 
one of them is chaotic, while the other is localized onto 
a classical torus. Thus the criteria for accepting that 
a quantum eigenstate corresponds to a quantized torus 
are the vicinity of the eigenenergy to the semiclassical 
prediction and localization of the eigenfunction. Based 
on these two criteria, we indeed identified 4 + 4 regular 
eigenstates (for the n = 0,1, m = 0...3 cases) in the 
spectrum. The accuracy of the adiabatic quantization is 
remarkable, since using the numerically-determined Ct 
for the semiclassical quantization shows that the differ- 
ence between the quantum and semiclassical energies is 
~ lO -1 ^^, except for the n = 0, m = 3 eigenstate [the 
leftmost one in Fig. 0a)], for which it is 0.835ns (here 
5ns = 5n/2 is the mean level spacing of the N-S system). 
If one approximates the curves Ct by semi-ellipses as in 
Eq. (0), the agreement remains the same for the states 
with quantum number m = 0, 1 and slightly deteriorates 
for the m = 2, 3 states. Finally, the observation that 



chaotic eigenstates are intermixed with regular ones in 
the given energy interval suggests that the Berry- Robnik 
conjectureSi for the spectrum of (normal) systems with 
mixed classical dynamics might also hold for ABs. 

Besides the regular eigenstates discussed so far, we 
have found that some of the bands of nearly constant 
T e h shown in Fig. 0b) with intermittent dynamics also 
support one or more quantum eigenstates. The energies 
of these eigenstates can also be obtained with semiclas- 
sical methods, although not as accurately as in the pre- 
vious case. We found that calculating the average time 
T e h of an electron-hole orbit in a given band and then 
using Eq. (0) with L = v^T e h one can usually predict 
the quantum eigenvalues with an error < 5ns- As an 
example, for the eigenstate shown in Fig. 0e) the quan- 
tum calculation gives e/A = 0.3202 while for n = 1 the 
semiclassical result is 0.3239, giving an error of 0.85ns- 
(There is also a quantum state which corresponds to the 
n = semiclassical one, but it is not as clearly localised 
as the one shown in Fig0e).) 

In summary, semiclassical analysis and the wave func- 
tion computation enable the classification of certain 
eigenstates as regular ones for Andreev billiards of mixed 
phase space. For regular states we present the first nu- 
merical calculation to show that EBK-like quantization 
scheme yields good agreement with the quantum results. 
Moreover, other states are either chaotic or can be asso- 
ciated with bands on the PS for which the time until the 
next Poincare section is approximately constant. 
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